function result=Fig1_f_2()
filename = 'gbsweep_277.GBs';
pathname = 'raw data/';

[data, ~] = load_GBs(filename, pathname);

n_bias = length(data.bias);
n_gate = length(data.gate);
peak = zeros(1,n_gate);
width = zeros(1,n_gate);
result.gate_2 = [];%pure
result.peak_2 = [];
result.gate_1 = [];%filtered
result.peak_1 = [];
result.valid = [];
width_2 = [];
result.nu = [];

for i=1:n_gate
    try
        tmp = fit(data.bias',data.map(:,i),'gauss1');
    catch
        tmp.b1 = 0; tmp.c1 = 1;
    end
    width(i) = tmp.c1;
    peak(i) = tmp.b1;
    
    result.gate_2 = [result.gate_2 data.gate(i)];
    result.peak_2 = [result.peak_2 peak(i)];
    
    
    if width(i)< 10e-3
        result.valid = [result.valid 1];
        result.gate_1 = [result.gate_1 data.gate(i)];
        result.nu = [result.nu 2.5 + (data.gate(i)-2.6525)/0.561 ];
        result.peak_1 = [result.peak_1 peak(i)];
    else
        result.valid = [result.valid 0];
    end
    
    
end

% % 
Fcolor='k';
Bcolor='w';
Fsize=18;
Bsize=14;
blur = 0.4;

f=figure();
plot(result.nu,-1*(1e3*(result.peak_1)+111.8),'LineWidth',2);

f.Color=Bcolor;
ax = gca;
ax.FontSize=Fsize;
ax.Color = 'none';
ax.XColor=Fcolor;
ax.YColor=Fcolor;
yl=ylabel('µ (meV)');
xl=xlabel('\nu');
yl.FontSize=Fsize;
xl.FontSize=Fsize;

x_limit=xlim;
% x_limit=[2.15,2.85];
y_limit=ylim;
y_limit=[-5,5];
axis([x_limit(1),x_limit(2),y_limit(1),y_limit(2)]);
end